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ABSTRACT 

We present an analysis of the luminosity distances of Type la Supernovae from the 
Sloan Digital Sky Survey-II (SDSS-II) Supernova Survey in conjunction with other 
intermediate rcdshift (z < 0.4) cosmological measurements including redshift-space 
distortions from the Two-degree Field Galaxy Redshift Survey (2dFGRS), the Inte- 
grated Sachs- Wolfe (ISW) effect seen by the SDSS, and the latest Baryon Acoustic 
Oscillation (BAO) distance scale from both the SDSS and 2dFGRS. We have analysed 
the SDSS-II SN data alone using a variety of "model-independent" methods and find 
evidence for an accelerating universe at >97% level from this single dataset. We find 
good agreement between the supernova and BAO distance measurements, both con- 
sistent with a A-dominated CDM cosmology, as demonstrated through an analysis of 
the distance duality relationship between the luminosity (d^) and angular diameter 
{(1a) distance measures. We then use these data to estimate w within this restricted 
rcdshift range (z < 0.4). Our most stringent result comes from the combination of 
all our intermediate-redshift data (SDSS-II SNe, BAO, ISW and redshift-space dis- 
tortions), giving w = — 0.81^Q']J(stat)±0.15(sys) and 17m = 0.22^Q°g assuming a flat 
universe. This value of w, and associated errors, only change slightly if curvature is al- 
lowed to vary, consistent with constraints from the Cosmic Microwave Background. We 
also consider more limited combinations of the geometrical (SN, BAO) and dynamical 
(ISW, redshift-space distortions) probes. 

Key words: (stars:) supernovae: general, (cosmology:) observations, distance-scale, 
cosmological parameters, large-scale structure of universe 



1 INTRODUCTION 

It is now widely believed that the late-time expansion 
of the universe is accelerating. General Relativity (GR) 
implies that the acceleration is driven by "dark energy" 
(DE) — an unknown energy component in the universe 
with a negative effective pressure. Describing dark energy 
by an equation-of-state parameter of w(z) — p/(pc 2 ) re- 
quires that w < —1/3. Alternatively, accelerated expan- 
sion could be an indication that GR is not the correct 
theory of gravity or that we have applied GR incorrectly 
in a cosmological context (see recent reviews of DE by 



Peebles fc Ratra 2003, Uzan 2006. JCopeland" et al hood , and 
Frieman. Turner, fc Hutererll2008l ) . 



Over the last decade, the most direct way of study- 
ing this acceleration of the expansion of the universe, and 
therefore DE, has been using Type la Supernovae (SNe), as 
they have been shown by many authors to be well-calibrated 
"standard candles" in the universe, i.e., their relative dis- 
tances can be determined from the dependence of their peak 
luminosity on the shape of the light curve. This method was 
used to great effect by astronomers in 19 98 to provide the 
first evidence for an ac cele rated universe l|Riess et al.l 1 19981 : 
iPerlmutter et al. I ll999l : see lFilippenkcll200a " for a review). 

Briefly, a Type la supernova occurs when a white dwarf 
star in a close binary system accretes enough mass from 
its compani on to undergo a ther monuclear exp l osion in the 
core. Both iPhillipj (|l993l ) and lHamuv et all (| 19931 ) have 
shown that such explosions can serve as consistent light 
sources in the universe to high accuracy. This is achieved by 
transforming the measured light curve of the explosion into 
the rest frame of the supernova (so called K-corrections) and 
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then correcting the luminosity at maximum as a function of 
the shape of the rest-frame light curve. 

Several techniques now exist for fitting S N light curves 
know n under different acronyms (Am is . lHamuv et al.l 
19961: MLCS. iRiess. Press, fc Kirshnerl Il996l: st retch^ 



Perlmutter et alj Il997l; CMAGIC, IWang etafl 120031 : 



BATM. iTonrv et all [2001 : SA LT, iGuv et all 120051 : AC12, 



Wang et all ^OOtF " SALT2, IGuv et all 120071 : SiFTO, 
Conlev et al1l2008l). In t hi s ana lys is, we consider MLCS2k2 
Jha. R icss. & Kirshnerl 120071 : |R iess. Press, fc Kirshnerl 



19961 ) . which is among the most commonly used, and best 



tested, of these various techniques. 

Recently, several dedicated SN surveys have been car- 
ried out to confirm and extend the earlier detect i ons of 
an acc elerating un i verse (HST, IRiess et alj |2004 l2007l; 
SNLS, lAstier et all 120061 : ESSENCE. IWood-Vasev et all 
120071) as well as new compilations of existing SN dataset s 
|Davis et ai]|2007l : iKowalski et~ai1l2008l : Ifficken et al]|2009h . 
In addition to supernovae, observations of Baryon Acoustic 
Oscillatio ns (BAO) can be used to measure distances in the 
universe dBlake fc Glazebrook 2003; S eo fc Eisensteinll2003l : 
iHu fc Haimanll2003l ). The BAO are caused by sound waves 
in the early universe which leave a preferred scale in the 
distribution of matter equal to the sound horizon at recom- 
bination. Today, this scale corresponds to ~ 100/ft Mpc 
(Hubble constant at present: Ho = lOO/i km/s/Mpc) and 
can thus be used as a "standard ruler" throughout the uni- 
verse. The BAO signat ure has been detected in the clus ter- 
ing of galaxy clusters bv lMiller. Nichol. fc Ba tuski (200l|), in 
the correlation of galaxi es i n the Sloan Digit a l Sky Survey 



ne cor relation 01 galaxies m 

SPSS lYork et allboOCh by [Eisenstein et alj (f2005h iHiitsil 
2006h . Ipadmanabhan et al. i fcoOTD . and lBlake et al.l(|2007 ' 



and i n the Two-degree F ield Galaxy Redshift Survey (2dF- 
GRS. IColless et allhoOll ) bv lCole et all (|2005l ). 



In addition to the geometrical methods discussed above, 
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observations of the dynamical properties of matter can pro- 
vide constraints on the matter density of the universe and 
w (assuming General Relativity is the appropriate theory 
of gravity). For example, the growth rate of structure in 
the universe can be observed via the coherent infall of 
galaxies into large clust ers and sup erclusters of galaxies 
seen in redshift surveys (|Kaiserl [l987T l. Also, the growth of 
structure can be measure d using the late-time Integrated 
Sachs- Wolfe (ISW) effect (|Sachs fe Wolfelll967D . which has 
now been detected to high significance through the cross- 
correlation of galaxy sur veys with the Cosmic M icrowave 
Background (CMB) (see iGiannantonio et al.ll2008l and ref- 
erences therein). The ISW is sensitive to deviations from 
a matter-dominated, Einstein-de Sitter universe (Qm = 1, 
where Qm is the matter density at present divided by the 
critical density.). 

Taken together, the present combination of cosmologi- 
cal measurements suggests we live in a flat univers, domi- 
nated by a cosmological constant (A) with the energy density 
in matter and A know n to a statistical acc uracy of better 
than a few percent fsee lDunklev et aT1l2009l ). However, sev- 
eral of these cosmological measurements, especially SNe, are 
becoming limited by t heir systematic uncer tainties which are 
now dominating, e.g., IHicken et ail (|2009l ) showed that the 
best combination of available SNe and BAO measurements 
provide 1 + w = 0.013lQQgg but with a systematic uncer- 
tainty of 0.11. Therefore, it is clear that future cosmological 
surveys must resolve these systematic errors through new 
observations and better analysis methods to mitigate their 
effect. 

This paper is one of three complementary papers fo- 
cused on the cosmological analysis of a new sample of inter- 
mediate supernova distances recently obtained by the SDSS- 
II Supernova Survey (see Section [5] for details). Our analy- 
sis differs from thos e pre sented in our companio n papers of 
iKessler et all (|2009h and lSollerman et all (|2009l ). as we first 
study the cosmological information obtained solely from the 
SDSS-II SN sample, and then in combination with other cos- 
mological p robes over the same redshift range (z < 0.4). Al- 
ternatively, IKessler et al. | (|2009h presents a detailed exami- 
nation of the impact of both statistical and systematic errors 
on deriving standard cosmological constraints based on the 
combination of the SDSS-II SN with most of the currently 
available high and low redshift SNe and which are all anal- 
ysed in a consistent way. ISollerman et al.1 (|2009h then use the 
same compilation of data to study an expanded set of exotic 
cosmological models, in combination with a wider variety of 
other cosmological information. Our approach is also com- 
plementary tothe many other analyses i n the literature (e.g . 
iDavis et al1l2007l . iKowalski et alj|2008l . IHicken et al.1 [20091 ') 
that have used data from all possible sources. 

In our approach we concentrate on the information from 
cosmological measurements that cover the same range of 
redshifts as the SDSS SN sample. Our aim is not to derive 
the most stringent limit on cosmological parameters avail- 
able, but rather to verify that if we restrict ourself to probes 
coming from a small and similar redshift slice the results 
on cosmological parameters remain stable and consistent. 
This approach is warranted because of the growing empha- 
sis on controlling systematic uncertainties in the analysis of 
cosmological data. There are clearly a number of system- 
atic uncertainties that could affect the use of SNe as cos- 



mological probes which likely depend on, or change with, 
redshift, including SN e volution (e.g., changes in the meta l 
licity of prog e nitor s t ars iTimmes. Brown, fc Truranl (|2003h 



iHowell et atl ilioo" 



dust (jConlev et al.1 \2QQ^ : 



Sullivan et al l d2009l )'l. intergalactic 
Holwerdal l 2008l ), Malmquist bias 



and t he effects of gravita tional lensing and peculiar veloc- 
ities (|Hui fc Greene! |200rj ). Moreover, the photometric un- 
certainties associated with combining SN data from multi- 
ple surveys, over a range of redshifts, is already seen as the 
main limitation in using presently available datasets (see 
IHicken et al.|[2009l ). Our analysis addresses this issue by fo- 
cusing exclusively on the SDSS SN sample, which is derived 
from a well-understood and stable photometric system. The 
SDSS has a relative photometric accuracy of better than 
2% in griz, and 4% for the u-band (Padmanabhan et al. 
2007), while the absolute calibration is also known to be of 
the order of 1%, leading to a homogeneous set of SN light- 
curve s with high photometric accuracy (see iHoltzman et alj 
I2008T ). This set of data is robust to uncertainties in light- 
curve fitting. For example, the MLCS2k2 and SALT2 fits to 
the S PSS-only sample ar e shown to agree well (see Section 
10 in IKessler et alTl2009h which is not the case when the 
higher redshift SN samples are added. 

The outline of this paper is as follows. In Section [5] we 
describe the SDSS-II SN data and use that data in Section 
3 to study the cosmic acceleration in the Universe. Section 4 
then compares the SDSS-II SN luminosity distances to the 
BAO distances from the SDSS and 2dFGRS, checking the 
distance duality relation. We then derive in Section 5 con- 
straints on w by combining the best-fit luminosity distances 
for SDSS-II SNe w ith the growth rate of structure measure- 
ments taken from lHawkins et alj (l2003f) and a new mea- 
surem ent of the ISW effect taken from IGiannantonio et all 
IMPOST ). We conclude in Section [6] 



2 THE FIRST- YEAR SDSS SN DATA 



The SDSS-II Supernova Survey (|Frieman et al.1 120081 ) was 
part of the SDSS-II project and was focused on construct- 
ing a large sample of intermediate-redshift SNe (0.045 < 
z < 0.42). One of the strengths of the SDSS-II SN Survey is 
that it builds upon the e xisting (and stable) i n frastructure 
from the original SPS S ( Fukugita et al.1 1 19961; Gunn et "all 



Eisenstein et al.1 



Strauss et al 



19981; I York : et al.lj2000l ; [Lupton et a~l200l(;lHogg et al.1 
200 if ; 



20021 ; IPier et al 



2001 



2003 



Gunn et al.ll2006l ). The SPSS-II SN Survey is based on re- 
peat imaging of "Stripe 82", a region of the original SPSS 
with significantly deeper photometry than the regular SPSS 
survey. This region is ~ 120° long and 2.5° wide, centered 
along the celestial equator and extending from 20 hrs to 4 rs 
in right ascension (passing through rs ). 

The SDSS-II SN Survey was carried out in three observ- 
ing campaigns from September through December in each of 
2005, 2006, and 2007 (there were also some observations for 
a short period in 2004) . The new imaging d ata were initially 



reduced using the standar d SDSS pipelines flS toughto n et al 
20021 ; ISmith et all |2002|; llvezic et al.1 [20041; iTucker et al 



2006 ; IPadmanabhan et al.lT2008l ; lAdelman-McCarthv et al 
20081 ). followed by specifi c image-subtract ion software to 
identify transient objects (|Sako et al.l|2008l ). To determine 
the nature of these transients, the data were both visu- 
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Figure 1. Residual Hubble diagram with respect to an empty 
universe for the 103 Type la SNe from the first year of operation 
of the SDSS-II SN Survey. The red line shows a ACDM model 
with (Oflf, ^a) = (0.3, 0.7), similar to our best-fit model given in 
Table [T] The blue line is the best fit to the data derived from the 
"sliding window" technique described in Section [4.1l and the cyan 
and green shaded regions correspond to the lcr and 2a confidence 
intervals, respectively. The black line indicates an expansion his- 
tory with a deceleration parameter of qo = —0.34 as described in 
Section [3] The lower panel shows the same parameterization but 
now converted to Dy according to Eq. (fOTt ■ The two data po ints 
represent the BAO measurements from iPercival et al.l J2009h . 



ally inspected and fit with SN models. Sub sequently, objects 
with a high probability of being a SN la (|Sako et al.1 120081 ) 
were spectroscopic ally observed using a variety of telescop es 
around the world (|Frieman et al.ll2008l ; Izheng et al.ll2008l ). 

In this paper, we only consider the 2005 observing cam- 
paign (the first year of operation) as the data from other 
years is still being collated and analysed. That year, the 
SDSS-II SN Survey discovered 130 spectroscopically con- 
firmed type la supernovae a nd an additional 16 spectro- 
scopically probable SNe la. In Kcssler et al.l (|2009l ) distance 
moduli are obtained for 103 of the spectroscopically con- 
firmed SNe la that pass stringent light-curve quality cuts, 
using the MLCS2k2 light-curve fitting routine. In the upper 
panel of Fig. [1] we show for these 103 SNe the residuals of 
the distance modules with respe ct to an empty unive rse. We 
refer the reader to Section 9 of iKessler et all (|2009T ) for an 
extensive discussion of systematic effects caused by changing 
the various light-curve fitting parameters. 

W e also direct the reader to Section 11 of lKessler et al.l 
20091) for a co mparison of the MLCS2k2 and SALT2 



Guv et al"1l2007l ) light-curve fitters. They show that for the 



SDSS-only data the systematic difference between these two 
light-curve fitting methods is only 0.04 in w, while Fig. 42 
of their paper shows the two methods give consistent dis- 
tance moduli for the same SDSS-II supernovae. We also 
highlight that t he two method s give s imilar contours in Fig. 
26a and 35a of IKessler et al.l (|2009l ) when comparing the 
full cosmological fits to the SDSS-only data. This motivates 
the analysis in this paper and means our results are unaf- 
fected by the choice of light-curve fitter used. We restrict 



th e analysis to the ML CS2k2 reduction taken from Table 10 
m IKessler et all (|2009l ), as these data include corrections for 
selection effects. 

When using these SDSS-II SNe for cosmological fitting, 
we calculate the confidence intervals via the x 2 statistic, 



2 

X = 



E 



(MLC,i(2») - Mmodcl(^,^)) 2 



+< 



(1) 



where and olc are the distance moduli and errors, re- 
spectively, derived from the light-curve fitting method (see 
IKessler et al"1l2009T ). and /imodei are the expected distance 
moduli according to parameters x of the assumed cosmo- 
logical model. Uncertainties in the measured spectroscopic 
redshift and peculiar velocity of the SN are taken into ac- 
count using 



0"sp 



(1 + *) 



ln(10) z(l + z/2) 



V(A 2 )2 + {Av p )*/<?, 



(2) 



where Az is the measurement uncertainty in redshift and 
Av p is the characteristic amplitude of the p eculiar veloci- 
ties, w hich we take to be 360 km s _1 (see, e.g.. lMasters et all 
l200rj ). We investigated the effect of correlated peculiar ve- 
locities of SDSS-II SNe and found no detectable effect (see 
Appendix B). We therefore do not include peculiar velocity 
correlations in further analyses. Besides peculiar velocities 
we ignore further potential correlations between individual 
SNe and treat them according Eq. {TJ as statistical indepen- 
dent. 

Following standard practice, we add an intrinsic disper- 
sion in distance modulus, of nt , in the denominator of Eq. (1) 
and determine it by setting xV-^dot = 1 for the best-fit cos- 
mological model. This term accounts for the fact that the 
errors on the distance moduls reported by the light curve 
fitter could underestimate the real error if we assume that a 
smooth cosmology is the correct underlying model. For the 
best-fit model with constant w, we fi nd (Tint = 0.088 mag for 
our SDSS-only SN sample, while in IKessler et al.l (|2009f ) a 
value of (Tint = 0.16 derived from the nearby supernova sam- 
ples is used in combination with the SDSS data. This results 
in broader contours compared to the one shown in this paper 
(see their Fig. 26a, and Appendix E for a possible explana- 
tion), but only marginally changes the most likely values of 
w and Q,m . Similarly if we omit any intrinsic dispersion we 
get narrower contours but only slight changes in w and Qm 
well within the errors given in Table [T] 

A further complication arises due to the uncertain cal- 
ibration of the absolute magnitude (at peak) of SNe la, 
leading to a degeneracy with the absolute value of Hq. To 
account for this, we marginalize over Ho- This procedure 
makes use of the relative distances reported by the light 
cur ve fitter but not o f their absolute value. In a recent pa- 
per iRiess et al.1 (|2009l ) redetermined the Hubble constant Ho 
or equivalently the absolute brightness of SN la. One could 
use this value obtained from measurements in the nearby 
universe as a prior on Ho in our analysis. But - as laid out 
in the introduction - we want to limit our analysis to probes 
taken at comparable redshifts as the SDSS SNe and there- 
fore refrain from using this additional information. 
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3 TESTING COSMIC ACCELERATION 

Given the homogeneity of the SDSS-II SN dataset, we begin 
our analysis by revisiting the original evidence for cosmic 
acceleration in the expansion of the universe. In detail, if we 
simply assume the universe is homogeneous and isotropic, 
and is described by a Robertson- Walker metric with a scale 
factor of a(t), then the purely kinematic deceleration param- 
eter, q, is defined by 



«(*) 



aa 

■ 9 ' 
a - 



(3) 



where q < corresponds to acceleration (we also assume 
light follows null geodesies and is therefore redshifted in the 
usual way). We relate q to the Hubble parameter by 



H(z) 



Hq exp 



i + g(fO 

1 + z' 



dz 



(4) 



Luminosity distance can then be calculated directly from 
the expansion history via 



d L {z) = c{l + z) 



dz' 



(5) 



where we have assumed a flat universe. Thus, the magni- 
tudes and redshifts of any SNe can be used to constrain q(z) 
without choosing a particular dark-energy model, or even a 
particular theory of gravity. The assumption of flatness is 
necessary in practice since constraints on q(z) degrade signif- 
icantly when curvature is allowed to vary. A prior on curva- 
ture from CMB measurements would, of course, strengthen 
the constraint, but such a prior is based on the validity of 
GR, counter to the intention of the q(z) analysis. 

The simplest deceleration model we can fit is a constant, 
q(z) = go- In this case, the luminosity distance simplifies to 



(6) 



Fitting Eq. JBJ to the SDSS-only SN data, we find a best 
fit of go = —0.34 and h = 0.636. Marginalizing the joint 
probability density function (PDF) over h, we find go = 
-0.34 ±0.18, or g < with 97% probability, i.e., the SDSS 
alone finds evidence for acceleration at 2a without concerns 
regarding the absolute calibration of the peak brightness of 
SNIa and the relative calibration between SN surveys. 

There is no reason to expect q(z) to be constant; indeed, 
under ACDM, q evolves from 1/2 during matter domination 
to —1 when vacuum energy dominates. However, we find 
that including additional parameters in q(z) (e.g. by Taylor 
expanding q) degrades our constraints to the point of being 
uninteresting. 

3.1 Principal Component Analysis 

When trying to reconstruct an unknown function from noisy 
data, there is the concern that particular features of the re- 
construction are not indicative of the true underlying func- 
tion, but an artifact of the chosen parameterization. This 
concern is magnified for a function like q(z), which is re- 
lated to the data fj,(z) by two derivatives. Therefore, we 
need a more robust way to determine if the universe has 
accelerated — a way that does not depend on the naive as- 
sumption of a constant q(z). A principal-component analysis 
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Figure 2. The principal components for q(z) derived from the 
SDSS-II SN data. Each mode includes a Fisher matrix estimation 
of the error bar, cr(ai), for its coefficient on. The errors are un- 
corrected, and we have ordered the modes according to the size 
of their error bars. Only the first six modes arc shown. 



( Hutcrcr fc Starkmanll2003l ) can be used to address this is- 
sue in a parameter-independent way. Principal components 
are a unique set of orthogonal basis functions (ei(z)), such 
that 



q(z) = auei(z), 



(7) 



which allows us to specify q(z) using the coefficients, o?i. 
The principal components, or "modes" ej(z), are explicitly 
constructed so that the coefficients can be measured inde- 
pendently of each other, i.e., they have uncorrelated error 
bars. To construct these modes, we start with a piecewise- 
constant parameterization of q(z) in bins of dz = 0.01, and 
we use our data to calculate a Fisher matrix for this param- 
eter set and the Hubble parameter Hq. After marginalizing 
over Ho, the modes ei(z) are given by the eigenvectors of the 
resulting matrix. We are free to normalize these functions 
so that 



ei(z)ej(z) dz = Sj, 



(8) 



which now specifies each function up to an overall sign con- 
vention. This procedure is completely general; in the limit of 
dz — > 0, we can specify any continuous function using these 
modes. The procedure is parameter-independent in the sense 
that we do not specify the modes a priori: they are deter- 
mined primarily by the data. We do choose q(z) = as 
the fiducial model for our Fisher matrix calculation, but the 
resulting modes are insensitive to this choice. 

The six modes constrained best by the SDSS-only SN 
data are shown in Fig. [2] Each mode is accompanied by an 
estimate of the error bar on its coefficient, cr(ai), calculated 
from the eigenvalues of our Fisher matrix. The first coef- 
ficient, ati, has the tightest constraint, while higher mode 
coefficients have increasingly larger errors and so provide lit- 
tle information about the shape of q(z). Note that the first 
mode does not cross zero, an d we have chosen its sign to be 
wholly positive. As argued bv lShapiro fc Turner! (|2006t ) . this 
mode is useful since in addition to being well constrained, 
Qi can only be negative if q(z) is negative for some z. There- 
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fore, if we fit the model q(z) = ctiei(z) to the data, and then 
measure ot\ < 0, it constitutes parameter-independent evi- 
dence that the universe has accelerated at some redshift. We 
need not consider additional modes since the constraint on 
qi is, by construction, independent of the constraints on the 
other at. Marginalizing over h, we find ai — —0.155 ±0.086 
and determine that ai < with 96% probability. The fact 
that the error bar from our fit closely matches the Fisher 
matrix estimate of a(cei) = 0.082 demonstrates consistency. 
This result is comparable to our qo fit above, while provid- 
ing robust evidence for an accelerating universe regardless 
of the model assumptions. 

The significance of our detection could be enhanced by 
combining the SDSS-1I SNe with other SN dataset s This 
has already been done in part by IShapiro fc Turner! (|2006h 
but could suffer from systematic uncertainties associated 
with combining data from different i nstruments and sur - 
veys. Again, the reader is referred to iKessler et alj (|2009T ) 
for a detailed description of such combinations of the SDSS- 
II SN Survey with other SN datasets. We plan to repeat this 
analysis with the full 3-year SDSS-II SN dataset. 



4 COMPARISON OF DISTANCES 

We next consider the comparison of our SDSS-II SN dis- 
tances with other geometrical distance estimates over the 
same r edshift range. This is motivated by the original re- 
sults of IPercival et al.l (|2007f ) who noticed some tension (at 
> 2a) between the cosmological constraints derived from 
nearby BAO measurements (z = 0.2 and z — 0.35) and 
higher redshift SNe of Astier et al. (2006). The BAO pro- 
vide a measure of distances in the universe by relating the 
scale of the sound horizon at last scattering (r s ) to the scale 
of the corresponding correlations seen in the galaxy distri- 
bution. One suc h measurement o f this ratio is given by the 
A-parameter in lEisenstein et al.l (|2005l ) . This parameter is 
frequently use d in combination with SN data to deriv e con- 
str aints on w (I Astier et al.|[200r3 ; iKowalski et alJl200Sl ) . e.g. 



see IKessler et all "! 2009h for the combination of the SDSS-II 
SN data with the A-parameter. 
Here, we adopt the quantity 



Dv(z) 



D M 



H(z) 



(9) 



defined in lEisenstein et al.l (|2005f) 
moving distance. IPercival et al.l 



where Dm is the co- 
[ 2007h showed by combin- 
ing measurements of r s /Dy at redshifts of 0.2 and 0.35 from 
both the 2dFGRS and SDSS galaxy samples, that one can 
obtain the ratio of the distance between two different red- 
shifts that is independent of both r 3 and Ho- This approach 
also avoids a large extrapol ation between the red shift of re- 
combination (zcmb = 1090. iKomatsu et alj|200"^ ) and these 
intermediate-redshift measurements. 

For the analysis presented in this paper, we have 
adopted the latest value of Dy(z = . 35) /-D y(z = 0.2) = 
1.736 ±0.065 taken from IPercival et al.1 (|2009h . The inferred 
value of D v (z = 0.35)/£V(0.2) from Percival et al. (2009) 
is lower than that of Percival et al. (2007), bringing it into 
better agreement with ACDM. This change was caused by 
a revised error analysis and a change in the methodology 



adopted, as well as the addition of more data. In this pa- 
per we do not use constraints on r 3 (zd)/Dy(z), which de- 
pend on the sound horizon at the baryon drag epoch r s {z&). 
We therefore avoid including CMB data, commonly used to 
model the sound horizon. 



4.1 The Sliding Window Method 

In Fig. [T] we show the Hubble diagram for the SDSS-II SN 
data discussed in Section 2 compared with a variety of cos- 
mological and non-parametric models discussed herein. We 
find a scatter of 0.14 mag around these models independent 
of the particular fitting method. The first non-parametric 
model we consider is a "sliding window" method which al- 
lows us to investigate the general shape and smoothness of 
the Hubble diagram without assuming a cosmological model. 
We have thus fit piecewise Hubble parameters and lumi- 
nosity distances in different redshift bin s using a local red- 
shift w indow following the approach of iDalv fc Diorgovskil 
(|2003h . At each redshift z, we fit the SNe co- moving dis- 
tances, Dm(zi) = di/fl + Zi), over a fitting window of 
Zi — Az < Zi < Zi + Az Q by a polynomial of second order 
given by: 



f(zi - z) = A + Ai(zi -z) + A 2 (zi - zf 



(10) 



The values of Ai are determined separately via x ~ 
minimisation in each redshift window, and we slide this win- 
dow as a function of redshift in increments of 0.01 through- 
out the entire range. The best-fit Dm at z is proportional to 
Ao, while the best-fit c/H is proportional to A\. and Dy is 
related to [zAoAi] 1//3 . Our results depend on the size of the 
redshift window, with a wider window allowing less flexibil- 
ity but smaller errors, and vice versa. 

We show in FigfT] the resulting non-parametric fit to 
the SDSS-II SN data as an example for a window size of 
Az = 0.15 which demonstrates that the SDSS SNe data 
is fully consistent with the individua l BAO measurements 
at z = 0.2 and 0.35 of IPercival et ah | (|2009h (the cyan and 
green shaded regions indicate the la and 2a errors, which 
are highly correlated as they share the same data points in 
the overlapping fitting windows). Next, we derive the ratio 
Dy(0.35)/ Dy (0.2) and determine the covariances between 
the Ai values within, and between, redshift bins using the 
observational errors. This is shown in Fig. for different 
window sizes. It is interesting to note that the sliding window 
method tends to prefer large values (steeper slopes) for Dy 
compared to Dy calculated from a qo fit described in Section 
3, or Dy from the best fitting w — Q. m parametrizations 
given Table 1. We also see that the sliding window method 
provides values for Dy which are fully consistent with the 
BAO result. We also show in Fig. [3] the best-fit ratios of 
Dy when the SDSS-II SN data are simultaneously fit with 
the BAO data. The SN constraints dominate these results 
because of their smaller uncertainty. 



1 We actually allow for a tapered region at both ends of this red- 
shift window where the weights (i.e., the inverse-squared errors) 
are reduced by a Gaussian function with a standard deviation 
of a z = 0.02. This leads to a suppression of fluctuations due 
to the inclusion of indiv idual data points into the window (see 
IDalv fc Diorgovskilliuol) . 
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Figure 3. Upper panel: Measurements of the ratio of 
Dy (0.35) /Dy (0.2) using different fits to the SN data, in com- 
parison to or in combination with the BAO (solid circle) mea- 
surement. The points labelled go indicate Dy derived using the 
fit in Section \3\ Points labelled w — Qm are derived from the 
u; — Qm parameterization marginalized over < Qjvf < 1 and 
—2 < w < 0. Points shown as triangles include the BAO measure- 
ment as a prior whereas boxes are without the BAO. The points 
to the right show the values derived from the "sliding window" 
method as a function of the window width. The grey points are 
for values of the sliding window that arc comparable in size to the 
redshift range of the whole SDSS data-set and therefore, approach 
the global fit. The dashed line indicates ACDM with Q.m = 0.3 
and = 0-7. Lower panel: Given the above-mentioned parame- 
tcrizations, we show the best-fit value of a defined in Eq. 11131 . 



4.2 Testing the Distance Duality Relation 

In the following we use the famous reciprocity relation 
|Etheringtonlll933l : lEllislll97ll ), or distance duality to com- 
pare the SN and BAO distance scales. In detail, the angular 
diameter distance and luminosity distance are related by 



£ = < I + *> 2 



("I 



(for a discussion, see e.g. iBassett fc Kunzll2004l ). This rela- 
tion relies on photon conservation, but holds for any geom- 
etry and any metric theory of gravity where photons follow 
null geodesies. Therefore, it is a general test of our underly- 
ing assumptions about the nature of our Universe. 

One might have expected that the distance duality 
relation has already been tightly constrained by observa- 
tions of the blackb ody CMB spectrum from the COBE FI- 
RAS experiment (|Mather et al.l Il994l ). However, this ob- 
servation does not constrain deviations from distance du- 
ality as the photon number may not be conserved (either 
through production or loss of photons) or more radically, 
photons may not follow null geodesies. Also, a grey dust 
component that absorbed photons independent of frequency 
would not cause spectral distortions away from a black- 
body in the CMB since all frequencies would be affected 
equally. However, this grey dust would cause strong devi- 



ations from distance duality since it would make the lu- 
minosity distance to any objects larger while leaving the 
angular diameter distance unchanged. Another way to hide 
the distance duality effects from CMB observations would 
be to affect photon number only at much higher or lower 
frequencies than the microwave. This is, for example, what 
was needed to make the axion-photon mixing proposal for 
the dimming of the SNe la consist ent with CMB constraints 
l|Csaki. Kaloper. fc Terning||2002l) . 

Beyond the CMB, several other analyses, using sim- 
ilar data to that discussed herein, have reported evi- 
dence for violations of distance duality at the ~ 2a level 
(IBassett fc Kunzl2004l : lLazkoz. Nesseris fc Perivolaropoulosl 
120081 ). We revisit this iss ue here using a methodology sim- 
ilar to that outlined by More. Bow, fc Hogg! (|2009l ) and 
lAvgoustidis. Verde, fc Jimenez! ( 20091 ) that does not rely on 
the absolute calibration of the distances to compute the 
ratio dh{z) / '(cIa(z)(1 + z) 2 ). Instead, we check the relative 
behaviour of this ratio as a function of redshift by testing 
the consistency of the ratio at two redshifts, z = 0.2 and 
z = 0.35 where we now h ave updated BAO measurements 
from iPercival et all (|2009T ). 

In the following we parameterize the distance duality 
relation in what we call the a-model as 



d L = (i + z y +2a d A = (i + z) 



(l+c) 



(12) 



where a = represents the expected distance duality rela- 
tion, and therefore a/(l indicates a possible violation. To 
quantify the discrepancy between the two measures, we re- 
place Dm in Eq. with dh from Eq. (I12|l . and derive the 
relation 



(l + z 2 ) I+2a zi _ D v ( Zl ) di(z 2 ) H(zi) 
(l + zi) 2 + 2a z 2 



D v (z 2 ) d 2 L ( Zl ) H(Z2)' 



(13) 



On the right-hand side, Dy (zi) / Dy (z 2 ) is given by the BAO 
measurements, while d 2 L {z 2 ) / d 2 L {z\) and H(zi)/H(z 2 ) have 
to be inferred from the SDSS-II SN data. We calculate dL 
and H(z) using the two methods introduced above. First, 
we use the results from the go fitting, and subsequently we 
calculate d^z, go) and H(z, go) at redshifts of 0.2 and 0.35. 
The parameter a, and its error, are then calculated by Eq. 
(|13|) . We note that Eq. [13] probes the consistency of the ratio 
at one redshift given the ratio at the other redshift and is 
therefore not sensitive to any scaling proportional to (1 + z) 2 
but would be sensitive to any other loss function. 

For the second method, we use the sliding window tech- 
nique to derive D v , S n = [z J 4o j 4i] 1/3 ( see Eq. JTD])) and 
the corresponding covariance matrix at the two redshifts 
(z = 0.2 and 2 = 0.35) where we have BAO measurements. 
The best fit and error of a is calculated by applying Bayes 
theorem. In detail we model Dy, bao = P{1 + z) 2a/3 Dy, 
at the two redshifts based on Eq. [13] where (5 is a free scale 
parameter absorbing Ho and the scale of the sound horizon 
r s at recombination. We then calculate the likelihood of the 
BAO Dy measurements for a and /3, by integrating over all 
possible SN Dy at z = 0.2 and z = 0.35 given the gaus- 
sian prior p(D'y\Dy : sn) constructed from the results of the 
sliding window technique: 



L(D V 



>\0,<*) 



dD'y X 
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Figure 4. Confidence intervals (68%, 95%, and 99%) in the w — 
!)m plane for a flat universe and a miCDM model derived from 
the SDSS SNe (shaded blue region), the BAO D v (0.35)/D v (0.2) 
ratio (filled red, upper 99% contours arc off the plot), redshift- 
space distortions (orange), and the ISW effect (green). In regions 
where the SN contours overlap the BAO contours the latter are 
indicated as red lines. 

[L(D v , BA oW(l + zf a/s D> v )p(D> v \D V}SN )] , (14) 

The covariance matrix used in L(Dv, BAo|y3(l + z) 2a ^ 3 D'v) 
is taken from Percival et al. (2009). Applying this procedure 
and subsequently marginalizing over p, we calculate the best 
fit a and its error. The results are shown in the lower panel 
of Fig. [3j along with the result from the go method. 

For the qo parameterization, we find a = —0.55 ± 0.45, 
while for the sliding window scheme, we find that all results 
agree with a = within one sigma, i.e., the errors on a 
are ~ 0.5 for most of the window sizes shown in Fig. [3] 
These results re-enforce our findings that the SN and BAO 
distance scales are now in good agre ement over this redsh ift 
range (confirming the new findings of lPercival et all (|2009h V 



5 CONSTRAINING COSMOLOGICAL 
PARAMETERS 

In contrast to the previous sections, which focused on kine- 
matic constraints of the cosmic expansion, here we inves- 
tigate the constraints on standard cosmological parameters 
using the SDSS SN data only in combination with dynamical 
measurements, eg. from the BAO, redshift space distortions 
and ISW. This will have less statistical pow er than using 
the co mbin ation of SN datasets p resented in iKessler et al.l 
(|2009h and lSollerman et all (|2009r i but our analysis is com- 
plementary to these companion papers and maximizes the 
impact of the homogeneity of the SDSS data (both for the 
SN and BAO). 

In Fig. [4] we begin by showing our constraints on w and 
Qm using only the SDSS SN data. This is achieved using the 
X 2 statistic according to Eq. {1} over a three dimensional 
grid of 200 bins in w, 100 bins in Qm and 160 bins in Ho 
ranging from 40 to 80 km s _1 Mpc _1 . Subsequently, we con- 
vert the x 2 into a likelihood using L = exp( — |(x 2 — Xmin))i 
where Xmin is the minimum x 2 found in the parameter space 
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Figure 5. Similar to Fig. [4] but now including curvature as de- 
scribed in Appendix IaI 

and, in our case, is by definition close to the number of SNe 
in the dataset (as we have added (Tint). We then marginal- 
ize over Ho by summing the likelihoods over the //o-bins. 
The shaded blue contours show the resulting confidence lev- 
els (68%, 95%, and 99%) in the w — SIm plane under the 
assumption of a flat universe (£Ia = 1 — Qm)- These con- 
tours include statistical errors only; systematic errors are 
discussed in Section f5. 2 1 

Fig. [5] is similar, but here - instead of assuming a flat 
universe - we allow for curvature according to the CMB shift 
parameter R. We do this by by calculating for a given (w, 
Qm) combination the corresponding value of 51a according 
to the constraints from R (see Appendix [X] for more details 
on the treatment of curvature). As discussed in Section 1, 
we see little effect on our results from allowing curvature to 
vary because of the relative small deviations from flatness 
allowed by the WMAP data. For results see Table 1. 

For comparison with the blue SN contours, we also pro- 
vide in Fig. [4] and Fig. [5] the (red ) cont ours for the BAO 
measurements from iPercival et alj (|2009h . The BAO mea- 
surements are in reasonable agreement with the SDSS-only 
S N contours but still p refer w < — 1 as originally discussed 
in lPercival et al.l (|2007l ). although all measurements are con- 
sistent with a cosmological constant. As the BAO and SN 
distances have reasonable overlap, we provide in Table[T]var- 
ious constraints on w and Qm derived from combinations 
of the SN data with intermediate redshift dynamical mea- 
surements described in the following Sections and the BAO. 
We note that the SN data dominates the width and posi- 
tion in the w-direction of the likelihood contours, as they 
have smaller errors compared to the BAO constraints, while 
adding the BAO data helps to curtail the large values of Qm 
seen for the SDSS-only results. 

5.1 Constraints from dynamical measurements 

To improve the cosmological constraints, we can include a 
number of other low redshift cosmological measurements 
thus providing a first measurement of the cosmological 
model within the local universe. In particular, we consider 
constraints on w derived from measurements of the growth 
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Figure 6. Confidence intervals (68%, 95%, and 99%) in the 
w — CIm plane for the combination of SDSS-II SNe with BAO 
(red), redshift-space distortions (orange), ISW effect (green) and 
for the combination of SN, BAO, RS and ISW (black) under the 
assumption of a flat universe. Numerical results are given in Table 



m 



of structure in the universe, including redshift-space distor- 
tions from the 2dFGRS and the ISW effect from the SDSS 
imaging survey. Both of these methods are particularly sen- 
sitive to Qm and thus provide independent and orthogonal 
constraints to the SN data. 



5.1.1 Redshift- Space Distortions 

Under the assumption that galaxies are related to the 
large scale dark-matter distribution, the anisotropy of the 
redshift-space correlation function depends on the parame- 
ter 



0(z) = f g (z)/b(z), 



(15) 



where b(z) is the linear bias relating the galaxies to the un- 
derlying dark matter and f g (z) the growth rate of structure. 
Absolute deviations between the real-space and redshift- 
space correlation functions depend on the parameter combi- 
nation f g (z)cr& (z), where erg (z) is defined as the root-mean 
square (rms) mass fluctuation in spheres of radius 8h _1 
Mpc, and provides a convenient way of no rmalising the mat- 
ter fl uctuations (for a recent review see iPercival fc White] 
I2008T I. 

To remove either the dependence on galaxy bias from 
a measurement of 0(z), or equivalently the dependence on 
o~s(z) from the measurement of f g (z)as(z), we need fur- 
ther cosmological information. In this paper, we adopt the 
central value of (3(z w 0.15) = 0.49 ± 0.09 calculated by 
lHawkins et~al1 (|2003l ) from th e 2dFGRS, wh ich is an update 
of the measurement given bv lPeacock et all d200lf) . We fol- 
low the procedure outlined bv lGuzzo et al.l (|2008h . and con- 
vert from /3 to f g by adding an additional uncertainty of 0.12 
in quadrature to account for the uncertainty in galaxy bias, 
which is estimated to be close to unity. This error includes 
the cosmological dependence of the bias measurement. For 
further calculations, we thus use f g (z = 0.15) = 0.49 ±0.15, 



Table 1. Results for w and Qm from combinations of the SDSS- 
II SN and BAO data with redshift-space distortions (RS) and the 
ISW effect. We provide measurements based on the MLCS2k2 
light-curve fitting method and our assumptions about curvature 
(Geo.). Values for w are derived after marginalizing over Qm and 
vice versa. The given uncertainties are statistical errors only; for 
an estimate of the systematic uncertainty see Section 15.21 



Dataset 



Geo. 



SN+BAO 


flat 


-0 74+ - 17 
u - '^-0.32 


u -°'-0.64 


SN+RS 


flat 


_ 77+0-19 

u - ' ' -0.25 


26+ 15 
U - ZO -0.12 


SN+ISW 


flat 


-0 74+ 16 

u - '^-0.22 


23 +0 - 1S 
U - ZJ -0.12 


SN+RS+ISW 


flat 


-0 76+ - 17 


23+° 10 
u - zo -0.08 


SN+RS+ISW+BAO 


flat 


(I O1+0-16 
-O-Ol-0 .18 


n 99 +0.09 
U - Z -0.08 



SN+BAO 


curved 


~0.99t° : » 


n 50+°- 13 

U.0U_ Q 22 


SN+RS 


curved 




31+ 015 


SN+ISW 


curved 


-0.781°^ 


27+ ' 15 
U - Z '-0.13 


SN+RS+ISW 


curved 


_Q 7Q+0-16 
u - ,y -0.20 


27+ ' 10 
u - z '-0.09 


SN+RS+ISW+BAO 


curved 


-0.85_ ; 19 


27+ ' 10 



assuming the weighted median redshift of the 2dFGRS sur- 
vey of z « 0.15. 

We adopt the parameterization given bv lLinderl (|2005l) . 



g(a) = — = exp / Gihia[S7ju(a) 7 — 1], (16) 
a Jo 

which is related to the growth rate by 

h = S ~, (17) 

where S — Spm/pm describes the perturbations in the den- 
sity of matter (p m ). For constant w > — 1, the exponent 7 
in Eq. (|16[> can be approximated as 



7 = 0.55 + 0.05(1 + w), 



(18) 



while for a phantom-like dark energy component (with w < 
— 1) the exponent is 



7 = 0.55 + 0.02(1 + w). 



(19) 



Solving Eq. (|17|l numerically, we derive the (orange) 
contours shown in Fig. [4] and Fig. [5] In Table [1] we present 
the constraints on w from the combinations of these data 
with the SDSS-II SN and BAO likelihoods, marginalized 
over Qm and vice versa. 



5.1.2 The Late-Time integrated Sachs-Wolfe Effect 

The integrated Sachs- Wolfe (ISW) effect is caused by the 
change in energy of CMB photons a s they pass through a 
time- varying gravitational potential (|Sachs &c Wolfej[l967l ). 
In a flat, matter-dominated universe, we would not expect 
to see an ISW effect as the large-scale gravitational poten- 
tials do not change in conformal time. However, in a uni- 
verse dominated by DE or curvature, we should detect a so 
called late-time ISW effect, which provides a direct measure 
of these quantities at the redshift of the changing potentials, 
i.e., the effect does not depend on the previous history of the 
growth of structure. 

The late-time ISW effect introduces additional sec- 
ondary anisotropies on top of the primary CMB fluctuations 
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and is therefore hard to detect directly. However, the ISW 
effect can be seen via the cross-correlation of the CMB 
with tracer s in the large-scale structu re of the universe as 
outlined by ICrittenden fc Tur ok (1996). This has now been 
achieved by a number of authors using a host of differ 



ent 
2002 


gaiaxy aarasets (rosaioa. uazranaea. &z oastanacr 
; Scranton et al.l 2003; iBoughn & Crittenden 


2004 


; Nolta et al. 20041; Afshordi. Loh. & Strauss! 


2004; 


Fosalba & Gaztanagal 2004; Padmanabhan et al.l 


2005; 


Cabre et al. 


2006; Giannantonio et al. 2006; Rassat et al. 


20071; Ho et all 20081; Granett, Nevrinck. & Szapudi 2008f). 



In this paper , we e xploit the recent analysis of 
iGiannantonio et al.l (|2008f ) and focus on the subset of 
intermediate-redshift (z < 0.4) SDSS data they used. Even 
this subset of data shows a dete ction of the ISW effect at 
the 3<t level l|Giannantoniol2008al ) . The contours for this new 
determination of the ISW effect are plotted in Figs. [4] [5] and 
[S] In Table [T] we present the combination of this new ISW 
effect measurement with our SDSS-II SN and BAO data, 
using the same procedure as discussed in Section [5.1.11 

5.2 Evaluation of SN Systematics 

We provide in Table [1] measurements of w and Qm 
from various combinations of the four data-sets consid- 
ered herein (SN, BAO, RS, ISW). The most stringent con- 
straint comes from the combination of all the probes giving 
w = -0.8llg:ig and Cl M = 0.22±° ™, which is competitive 
given the restricted redshift range considered in this analy- 
sis. However, much of this constraint comes from the com- 
bination of just the SDSS-II SNe and ISW measurements 
(See Table [1]). The ISW contours already include correla- 
tions between different angular and redshift bins and cosmic 
variance, and could there fore be considered stable (see e.g. 
IGiannantonio et al.ll2008T ). Similarly the contours we use for 
redshift-space distortions include the dominant uncertainty 
coming from t he galaxy bias (see t he procedure laid out in 
Section T5. 1.1 p . IPercival et al.l (|2009l ) has done several checks 
and found that their result is robust against variations in 
sample selection, number of redshift slices, calibration and 
other potential influences. Therefore, the results presented 
in Table 1 includes major uncertainties affecting the other 
probes but only the statistical uncertainties from the SDSS- 
SN data on the mea sured cosmological p arameters. 

As discussed in iKessler et~alT(|2009l ) . the SDSS-II SN 
distances also depend on the detailed choices and assump- 
tions within the MLCS2k2 supernova light-curve fitting pro- 
cedure, including different training vectors, priors on Av and 
Rv, uncertainties in zero points and the filter systems, and 
selection biases. To quantify the systematic uncertainties as- 
sociated with these parameter choices, we repeat our analy- 
sis above for th ese different choices and, following the proce- 
dure laid out in IKessler et ail (|2009h , we calculate a variation 
of Aw = ±0.15 with respect to the fiducial MLCS2k2 reduc- 
tion presented in Table [1] in the case of combining all four 
constraints and slightly larger values for the other cases. 

Our estimates of the systematic uncertainty for 
MLCS2k2 a r e lar ger than the values calculated in 
IKessler et all (l2009h bec ause t hey use the BAO ^-parameter 
from lEisenstein et al.l (|2005T l, and the added constraints 
on Qm and w derive d from using the CMB i?-parameter 
jKomatsu et al.l I2008T ) . We reproduce their values for the 



systematic uncertainties on MLCS2k2 (Aw w 0.1), if we in- 
clude these constraints in our analysis. However, in this pa- 
per, we restrict our analysis to intermediate-redshift probes 
and therefore do not include the CMB constraints, which 
results in larger uncertainties. 

Our analysis of the MLCS2k2 systematic uncertainties 
discussed above does not inc lude the large shift in w dis- 
cussed in IKessler et ail (|2009l ) when the rest- frame [/-band 
template is removed in the l ight-curve fitting. As seen in 
Table 6 o f IKessler et all i|2009l ) . removing the rest-frame U- 
band results in a —0.31 shift in w, while we find a shift of 
—0.43 if we remove this data. This particular uncertainty 
would therefore give rise to a bimodal result either centered 
around w ~ —0.8 (with U-band included) or w — —1.2 
(without U-band), yet both consistent with w = — 1 within 
the errors. 

We do not add the uncertainty due to excluding the 
rest-frame U-band to our systematic errors because we be- 
lieve it is incorrect to exclude this data from the light -curve 
fitting. Even though there is evidence for diversity in the UV 
spectra of SNe la (see Ellis et al. 2008; Foley et al. 2008), the 
removal of the rest-frame U-band data from the SDSS-only 
analysis results in the light curve fitter using only two filters 
at z < 0.2 to constrain the colors of the SNe. This provides 
significant freedom to the MLCS2k2 fitter and then the pri- 
ors on the fitted parameters become important. We note 
that w is only shifted by —0.1 when using t he SALT light 
curve fitter (see Table 8 in lKessler et~aT1l2009l ') on the SDSS- 
only sample with the rest-frame U-band excluded. This is 
the only noticeable difference between these two light curve 
fitting methodologies when considering the SDSS-only sam- 
ple; namely the e rror on w wh e n the rest-frame U-band is 
excluded. Finally, IKessler et ail (|2009l ) also sees a clear jump 
in the SDSS Hubble diagram at z ~ 0.2 when the rest-frame 
U-band is excluded from the MLCS2k2 analysis (see their 
Section 10.1.3 and Fig. 30), indicating that a constant w 
model is not a good fit in this case. 



6 CONCLUSIONS 

We present an analysis of the luminosity distances of Type 
la Supernovae from the Sloan Digital Sky Survey-II (SDSS- 
II) Supernova Survey in conjunction with other intermedi- 
ate redshift (z < 0.4) cosmological measurements including 
redshift-space distortions from the 2dFGRS, the ISW ef- 
fect, and the BAO distance scale from both the SDSS and 
2dFGRS. We have analyzed the SDSS-II SN luminosity dis- 
tances using several 'model-independent' methods, includ- 
ing fitting the data using a q(z) parameterization, prin- 
cipal components, and a non-parametric "sliding window" 
method. We find consistent results between all these meth- 
ods that provides evidence for an accelerating universe based 
solely on the first-year SDSS-II SN data. The strongest ev- 
idence we find comes when we make the strongest assump- 
tions, that go is constant and the universe is flat which gives 
probability for acceleration of > 97%. 

We also compare our SDSS-II SN data with the local 
BAO measurements, and find they ar e in good agreement . 
This is in contrast with the findings of lPercival et al.l (|2007f ) 
which found tension between the two distance measures, but 
confirms the new BAO analysis of Percival et al. (2009) who 
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note that this tension has now lessened. Taking this obser- 
vation further, we test the distance duality relation, i.e., for 
any metric theory of gravity, we expect di/(dA(l±z) 2 ) = 1. 
We see no evidence for a discrepancy from this relation (at 
the one sigma level) in contrast to previou s claims for a po- 
tenti a l violation on the 2a level as seen in dBassett fc Kunzl 
12004 : lLazkoz. Nesseris fc Perivolaropoulosl |200S|) . Finally, 
we present a new measurement of the equation-of-state pa- 
rameter of dark energy using a combination of geometrical 
distances in the universe and estimates for the growth rate of 
structure. Our strongest constraint comes from the combina- 
tion of all four data-sets discussed herein (SDSS-II SN, BAO, 
redshift-space distortions, ISW) with w = — 0.811q '^g (stat) 
and Qm = 0.22^®'^ (stat) (assuming a flat universe). How- 
ever, the combination of just the SDSS-II SNe and the ISW 
measurements alone is almost as powerful in constraining 
these parameters (Table 1). Our results only change slightly 
if we allow curvature to vary, consistent with the CMB mea- 
surements (see Appendix We quote a systematic uncer- 
tainty of Aw = ±0,15 based on the detai ls of the MLCS2k2 
light-curve fitter fsee lKessler et aL I l2009l for a fuller discus- 
sion) . 

Thus we have shown that low-redshift cosmological 
probes give a self-consistent picture of the distance-redshift 
relation. When combined with growth of structure and ISW 
at the same epoch that picture is consistent with ACDM 
and re-enforces the complementarity amongst other data 
and analyses in the literature. 
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Figure Al. The 68% contours in the £7 a — plane derived 
from the "scaled distance to recombination" R (see text) taken 
from the WMAP-5 results, or R = 1.715 ±0.021. The dot-dashed 
line indicates a fiat universe. 
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APPENDIX A: COSMIC CURVATURE 

In most parts of this paper, we assume a flat universe, or 
Qm + Q\ = 1. Observationally, the most robust constraints 
on curvature come from the distance to the last scattering 
surface of the CMB as determined by measurements of the 
CMB power spectrum. In cases where we consider curvature, 
allowing any combination of Qm and Qa , we include a prior 
on the "scaled distance to recombination" R defined as 



R= \J Qm H{ 



rfi(-ZCMB) 



(Al) 



°(1 + ^cmb)' 

given m IWang fe Mukheried with a value of R = 

1.71 ± 0.019 from lKomatsu et al.l (|2008l ). We note that R is 
independent of Ho because <1l scales linearly with 1/Ho. For 
fixed values of Qm and w it is now possible to constrain Qa 
from a measurement of R. We show 68% confidence levels in 
Fig- ED calculated from Eq. (|A1[) . As R does not depend on 
the Hubble constant Ho, or on the baryon density Qbh 2 , no 
further assumptions are needed on these quantities. From 
Fig. IA1I it is obvious that using only R, a curvature of order 
a few percent cannot be neglected (depending on the value 
of w). We introduce the effects of curvature in the redshift- 
space distortion analysis by adding an additional \ 2 term 
in the likelihood analysis calculated from Eq. (|A1[I and sub- 
sequently marginalize over Qa, where as for the ISW effect 
we pick for a given (w,Qm) combination the best-fit value 
of Qa- 



leading to significant covariance between pairs of SNe, i.e., a 
pair of SNe at radial positions , r,- have a projected velocity 
correlation function of r.,) = ((v(rj) •fj)(v(rj) •? j)). We 
can calculate this function in linear theory using the matter 
power spectrum, the linear growth function and its deriva- 
tive. This is interesting because, if this effect is detected, it 
may enable the SNe to constrain the parameters of structure 
formation, in addition to the standard background expan- 
sion. 

The e xpression for the full covarian c e between SNe 
is give n by iGordon. Land, fc Slosarl (|2007l ) ; lAbate fc Lahavl 
|2008l ) 



+ 

Hdr 



(B2) 



This can be compared with the standard diagonal ran- 
dom errors, which are 

,..2 /lnl0\ 2 r2 ^ /^i.A + 

where the intrinsic magnitude and velocity scatters a m ,a v 
have been introduced as usual. A numerical evaluation shows 
that the two are comparable at low redshift. In particular, 
for a pair of SNe at z — 0.05 and zero angular separation, 
the covariance is C v (i,j) ~ 0.1a(i)a(j). This decreases at 
higher redshifts and greater separations. 

This effect has been detected by 

IGordon, Land, fc Slosarl (l2007tl using a catalo gue of 124 low 
redshift SNe bv lJha, Riess. fc Kirshneil (|2007t l at z = 0.017, 
and it has been carried further t o constrain paramete rs 
such as as and the growth factor 7 (| Abate fc Lahavll2008l ). 

Here we repeated the analysis for the SDSS SNe, but 
since our minimun redshift is z ~ 0.05, we expect the ef- 
fect of correlated PVs to be small. Indeed, we found that a 
likelihood study performed by a Monte Carlo Markov chain 
(MCMC) analysis of the cosmological parameters yields no 
change in the results whether the full PV covariance matrix 
of Eq. (|B1[) is included or not. For example, if we set the 
intrinsic scatters a m — a v = 0, we find that the reduced 
X 2 jv for the best fit cosmology decreases by only 1% when 
the PV covariance matrix was included. 

Therefore, we are unable to detect the correlation of SN 
peculiar velocities with this data and are safe to ignore them 
in further analyses. This effect will become more important 
with larger samples of low redshift SNe. 



APPENDIX B: CORRELATED PECULIAR 
VELOCITIES 

The peculiar velocities (PVs) of supernovae (SNe) introduce 
an additional scatter onto t he Hubble diagram (s ee Eq. @). 
However, as pointed out bv lHui fc Greene! (|2006l ). we expect 
these PVs to be correlated, especially at low redshift, thus 



